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Abstract 

It was recently suggested that the gravitational action could contain a scale-dependent 
' cosmological term, depending on the length or momentum scale characteristic of the 
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processes under consideration. In this work we explore a simple possible consequence 
of this assumption. We compute the field generated in empty space by a static spherical 
source (the Schwarzschild metric), using the modified action. The resulting static potential 
turns out to contain a tiny non-Newtonian component which depends on the size of the 
test particles. The possible relevance of this small correction for the analysis of the recent 
Pioneers data [J.D. Anderson et al., Phys. Rev. Lett. 81 (1998) 2858] is briefly discussed. 

04.20.-q Classical general relativity. 
04.60.-m Quantum gravity. 
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1 Introduction 



On the base of the work by Hamber and Williams in Euclidean lattice quantum gravity |T|, it was 
recently suggested that the gravitational action could contain a scale-dependent cosmological 
term . In terms of quantum field theory this means that the effective value of the cosmological 
constant would depend on the characteristic length or momentum scale of the processes under 
consideration. Obviously since quantum gravity is not renormalizable, at least in the pure 
Einstein formulation, the concept of an effective coupling depending on the scale is not well 
defined at the perturbative level. A non perturbative analysis is needed, like that of Ref. 

The issue of a scale-dependent cosmological term is connected to the well known puzzle 
of the "global" cosmological constant A. In principle this can be present in the gravitational 
action already at the classical level and is expected to receive a large contribution from the 
vacuum fluctuations of the quantum fields; still, its observed value is very close to zero. Several 
theoretical conjectures were formulated in order to explain this cancellation ||. The experi- 
mental upper limits on a global A 0] stem from observations at cosmological scale and also 
from astronomical measurements in the solar system, where a negative A would correspond to 
a - never observed - Yukawa range for the gravitational potential of the sun and the planets. 

1.1 The scaling law 

How does the effective value A e ff of the cosmological constant (i.e., the average value of the 
scalar curvature evaluated in a certain 4D region) depend on the scale in the quantum lattice 
theory? It turns out that A e ff tends to zero when the size £ of the region goes to infinity, but 
for any finite value of £ it is given by the power law 



where I is the lattice spacing and 7 is a critical exponent. This scale dependence makes sense 
as long as we regard / as a small but finite quantity. In fact in lattice quantum gravity, unlike 
in the usual lattice theories, the lattice spacing not only acts as a regulator, but also represents 
the minimum physical distance. Thus I cannot be smaller than the Planck length and in the 
following I will indeed denote lp ~ 10~ 33 cm. (In this work we employ natural units, such that 



h — c — 1. Since in terms of h, c and G one has lp ~ JGh/c 3 , the scale dependence of the 



cosmological term can be regarded as a quantum effect.) 

Note that being I extremely small, the ratio //£ is also very small at any reasonable scale. 



(|A e// |G) (fl~(i/0 
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Also note that in order to isolate |A e jj| in eq. (|1|) we must do some hypothesis on the behavior 
of the effective Newton constant G. It is generally believed - and the numerical simulations 
confirm this to some extent - that the scale dependence of G is quite weak. Thus after extracting 
|A e //| from the adimensional product |A e //|G one finds that |A e jj| varies essentially as 

The critical exponent 7 has been numerically estimated in the simulations, but only for 
small lattices. Some consistency arguments (compare Section 2) allow to conclude that 7 is 
slightly larger than 2. The simulations also show that the sign of A e ff is negative: when the size 
of the 4D region, where the average of the effective value of A is computed, tends to infinity, 
A e ff tends to zero from below. 

This feature is familiar in numerical quantum gravity on the Regge lattice, where a positive 
value of the effective A always corresponds to the unstable, collapsed phase of the system. It is 
also connected, in the weak field limit, to a property of the linearized theory, where a positive 
cosmological term corresponds to an imaginary mass for the graviton and thus to unstable 
solutions. 

Although most of the considerations above refer (like the numerical simulations) to the 
Euclidean field theory, the assumption is made that in the weak field limit this is equivalent 
to the theory in Minkowski space. However, even if we do not trust this equivalence, we can 
regard eq. ([[]) as a reasonable ansatz for the scale dependence.^ 

1.2 Physical implications 

Suppose we accept the possibility of a scale-dependent cosmological term, stronger at short 

distances than at large distances. The next step is to understand the physical implications of 

this formal assumption. Which physical processes could be affected by such a term? How can 

it be taken into account? 

In an earlier work we studied whether a local cosmological term could let a graviton decay 

into collinear gravitons of smaller frequencies. The decay of a massless particle is a quite exotic 

process, whose amplitude obeys certain general rules |J; one of these is that decay is possible 

if the massless field self-interacts with a coupling of positive mass dimension (like A for the 

gravitons). However, for any field model comprising a scale-dependent coupling, the separation 

2 Bcsides the numerical simulations, there is a further argument showing that a small negative local cosmo- 
logical term with the scaling law (|l|) is necessary for the stability of the gravitational vacuum: it stabilizes the 
Einstein action with respect to a special set of field configurations, called "zero- modes" |5|. This holds in the 
Euclidean as well as in the Minkowskian formalism. 
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of long distance and short distance effects is not trivial. In the case of gravity this problem 
is even more serious, because while speaking of gravitons as of elementary particles we rely 
on the concept of Lorentz symmetry; but at the same time we admit that a decay of these 
particles could be caused by a non vanishing cosmological term, which corresponds classically 
to a curved background. Therefore it was suggested to consider separately (i) the phase space 
for the decay products at large distance (referred to flat spacetime), and (ii) a local decay 
amplitude (where the effective value of A would appear). This program, however, faces too 
many technical difficulties and inconsistencies to allow a clear resolution. 

This work aims at investigating the possible effects of a scale-dependent cosmological term 
in a simpler, classical context. In Sections 3 and 4 we study the motion of a body with given 
finite size in a background gravitational field, under the assumption that this body also "feels" 
a nonzero local cosmological constant, depending on the body size as in eq. ([!]). We wonder 
whether the motion is affected, and to what extent. 

There are some paradoxes and physical mistakes which lie in wait for us along this calcu- 
lation. The first, serious problem is the violation of the equivalence principle. In the presence 
of purely gravitational forces, any dependence of the center of mass motion of a test body on 
its size represents such a violation and could be in principle detected by an Eotvos-like experi- 
ment .0 Decisive points are, of course, the magnitude of the violation and its exact dependence 
on the test body size. This dependence could be so weak to lead to negligible differences in the 
motion, except for bodies differing in size by several magnitude orders.0 

A further crucial problem is a possible violation of the Lorentz symmetry of free motion. 
In an empty space, how could the motion of a test body be affected by the local A-term 
without spoiling Lorentz symmetry? Certainly the body cannot "slow down" or "turn" in 
some direction. Our calculations show that Lorentz symmetry is actually preserved and that 
the A-term only has an effect in the presence of a non-flat background defining a preferred 
direction. 

In order to fix the ideas we chose a Schwarzschild background metric, i.e., the static spher- 

3 According to eq. (|l|) and the following discussion, the scale dependence of A e // is a quantum effect, an 

"imprint" of the Planck length physics. One could dispute, at this point, how the quantization of a theory 

intrinsically invariant under diffcomorphisms can lead to results in contrast with the equivalence principle. In 

fact, as already discussed by Hamber and Williams, the quantum version of General Relativity on the Regge 

lattice reproduces the complete diffcomorphisms invariance only in the long distance limit. 

4 The experimental data of Ref. [Q are quite inspiring under this respect. Compare our discussion in Section 

5. 
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ically symmetric field generated by a certain mass M. We found an expression for the modifi- 
cation of this metric due to the inclusion of local cosmological term in the gravitational action, 
and studied the motion of a test body in the resulting field. We concluded that the main effect 
of the cosmological term is just a small unobservable "renormalization" of the product GM, 
independent on the test body size; in addition there is, however, a tiny non-Newtonian force, 
larger for smaller bodies and depending on the velocity v(r) (but unrelated to the familiar 
post-Newtonian corrections) . 

The discussion of this puzzling result is the subject of Section 5. 



Before embarking in the main calculation, let us derive in this section an estimate of the 
exponent 7 in eq. (p]). This is obtained just imposing a phenomenological consistency condition. 
From eq. (JTJ) we can isolate |A e //|(£), obtaining 



The range p of a Yukawa component in the gravitational potential is proportional to \A e ff 
We thus have (up to a numerical factor of order 1 which is already present in eq. (0)) 



We see from this equation that if 7 was exactly equal to 2, then p = £, i.e. the range of the 
Yukawa component of the gravitational field observed at the scale £, would be equal to £. 
But we know that such a Yukawa component is never observed, thus we must have p ^> £. 
For instance, by observing the planetary motion we can set an upper bound on the average 
curvature of spacetime at a scale of the order of the solar system size. 

We therefore assume that 7 is slightly larger than 2 and define a = 7/2 — 1, where a is a 
small positive number (of the order of 0.03 or less - see below). In this way, eq. @ becomes 



2 The critical exponent in the A e //(£) law 



|Ae//|(£) ~r 2 r 7 F. 



(2) 



p „ jl-7/2£7/2_ 



(3) 



P 




(4) 



and the ratio between p and £ is given by 




(5) 
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This ratio does not actually need to be very large, because the ratio between a possible Yukawa 
term Vy in the gravitational potential and the usual Newton term V N reduces to an exponential: 

7r~ ex p(-)- ( 6 ) 
Vn \pJ 

Thus p/£ is the logarithm of the ratio between the upper bound on Vy and Vat. Taking, for 
instance, Vy/V^ = 10~ 9 , one finds p/£ = ln(10 9 ) ~ 21. The size £ of the solar system is of 
the order of £ = 10 10 Km = 10 15 cm, thus the ratio £/Hn eq. (||) is of the order of 10 48 . By 
requiring (10 48 ) a < 21, one finds a < 0.03 and 2 < 7 < 2.03. 

3 The correction to the Schwarzschild metric 

In this section we compute to lowest order in the weak field approximation the field generated 
by a pointlike static source, assuming that the field equations contain a cosmological term 
localized at a given point. 

3.1 Einstein field equations and harmonic gauge propagator 

Let us start by fixing our conventions. The Einstein equations in the absence of a cosmological 
term are 

R^u{x) - ^g^{x)R(x) = S-kGT^x). (7) 

The trace of this equation is R(x) = 8tcGT(x), where T is the trace of T^ u . 

We focus on the case of a weak field g^ix) = r/^ + h^x); neglecting terms quadratic in 
h, one obtains the usual linearized version of Einstein equations, which in harmonic gauge can 
be written in the form 

K^ pa {x)h prj {x) = T„ v {x) (8) 
where K is a linear differential operator. Given the source T^, the solution of this equation is 



M*) = J ^vP^{x,y)T^(y) (9) 
where P is the propagator, that is, the inverse operator of K, given by 

W*.y) = -— ? Q T ( 10 ) 



with Q a constant tensor: 



In all these equations it is understood that, in accordance with the linearized approximation, 
the indices are lowered and raised with the metric 7]^ u , which is diagonal and has signature (- 
1,1,1,1). The square of a 4-vector x is x 2 = — (x ) 2 + x 2 . 

3.2 Field equations with a pointlike source and a local cosmological 
term 

In the presence of a cosmological constant A independent on x, eq. (0) takes the form 

RfMu(x) - ^g fll/ (x)R(x) + ^Ag fiu (x) = -8nGT^ u (x). (12) 

If = 0, this equation admits a solution with constant scalar curvature R(x) = 2A. 
The linearized version of eq. ( |12"D is 



- ^GK^ p(J (x)h»°(x) + [ Vfiu + h^(x)} = -SttGT^x). (13) 

We want now to solve this equation in the case when the energy-momentum tensor is that 
of a pointlike mass placed at the origin of the coordinates, called T^ u (x), and the effective 
cosmological constant depends on x (A — > A e ff(x)). Let us decompose the metric h pu (x) as 
follows: 

VW = h%{x) + A eff (x)h^(x) (14) 
where the metric h s satisfies the equation 

K^ pa {x)h s ^{x) = 7£(x) (15) 

and can be therefore regarded as the linearized version, in harmonic gauge, of the Schwarzschild 
metric. (We shall derive h s in a moment.) By introducing the decomposition ([14]) into eq. ( |i~3|) 
we obtain 



torGK^x) \h s ^{x)+A eff {x)h A ^{x) 



+ 



1 



+ 2 A efA X ) + k M X ) + A efA X ) h U X )\ = -^GT^(X). (16) 

Now let us write the function A e //(x) in the form A e ff(x) = A e fff(x), where f(x) is an 
adimensional function vanishing outside a region of width d (size of the test particle) centered 
at X (center of mass coordinate of the test particle), and such that /(X) = 1. In eq. flltf ) we 
can single out terms of zeroth order in A e fj, which cancel out due to eq. (|15|), one term of order 
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A^j, which can be neglected, and terms of the first order in A e ff, which lead to the following 



equation for the metric h A : 



1 



- 8nGK, upa (x)f(x)h A ^(x) + -f(x) [v,u + h%(x) 



0. 



(17) 



If we know the background metric h (x), we can find the correction h A (x) due to the 
cosmological term in the neighbourhood of X by applying the operator P, inverse of K: 

KM = 8n G f{x) J d\P»v pa {x,y) l -f{y) \f + h s ^{y) 



3.3 Computation of h s 

In order to find h s (x), we replace T (x) by its explicit expression 

T s (x) = M5 3 (x)5 p0 8 u0 

and integrate eq. (§): 



KM 



d 4 yP flupa (x,y)M6 3 (y)6^6 

M J dy°P tlvOO (x;y ,y = 0) = 
-2G 



M- 
2GM 



(2?Vo^o - VnvVoo) I dy 



(x° - y ) 2 - xe 



Then, in conclusion, h s is given by 



^oo (x) 



2GM 



2GM 



'no J 



(19) 



(20) 



(21) 



|x| |x| 

which is correct, since in general in a static field one has h 00 = —2V Newtonian \ the ii components 
are the same: 



ASM 



2GM , . 2GM 

—r-T {-VaVoo) = —r-T- 



(22) 



3.4 Computation of h A 

Now we can compute h A . In particular, we are interested into the 00 component, which gives 
the correction to the static potential felt by the test particle. We evaluate it at the center of 
mass coordinate X = (X Q , X) and find 

Ko(X) 

1 r ,a 2r] 0p r] 0a - T] 00 r] pa 



dry- 



(X - yf - is 



f(y) V ^ + h s ^(y) 



(23) 



Before looking at the integral, let us focus on the algebraic part of this expression. The term 
rf u in the bracket r\ pu + h s ' pa (y) , after multiplication by (2r] 0p r] 0a — r] 00 r] pa ), gives 

(2r] 0p r] 0a - r] 00 r] pa ) rf 7 = 2 Voo - 4 Voo = -2 + 4 = 2. (24) 

In a similar way, for the term h s ' pa (y) we obtain 

= 2/4(y) - (voo) 2 h s 00 (y) - mo E h Uv) = 

% 

2GM . x 8GM , , 

= —— 2-1 + 3 ) = - i - r . 25 

|y| |y| 

Now let us write the function / in the following form 

f( y ) = 1 if | y - X| < d; f(y) = elsewhere (26) 

or 

f(y) = e(d-\y-X\). (27) 
Setting for simplicity X = 0, since everything is static, we obtain the correction 
,,A m _ 1 />„ »(d-\y-X\) ( 2 SGM\ 

Introducing the new integration variable z = X — y, we obtain 

; '°°< x > = - a? / * / fe z* - - 5 ( 2 + N^q j ■ (29) 

By integrating the first term in the bracket (2 + 8GM/|z — X|) one clearly obtains a constant 
independent from X; this constant is irrelevant in the present case, since it does not give any 
contribution to the force. (This was expected: in flat space there cannot be any force on the 
test particle due to the local cosmological term, otherwise Lorentz invariance would be spoiled.) 
From the second term in the bracket (2 + 8GM/|z — X|) we obtain 

a , v , GM r r , 6(d- Izl) 1 

Since 

/\ IT 
z 2 — (y ) 2 — ie |z| ^ 

(compare the calculation of /i 5 ), we finally find 

C( X)^/, Z ^=M. ,32) 

7T J Z Z — A 
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This integral is well known: apart from the factor — 1/ir, it is the Newtonian potential generated 
at the point X by a mass distribution with spherical symmetry centered at the origin, having 
density p(|z|) = 6(d — |z|)/|z|. The integral is convergent, positive, and proportional to 1/|X|; 
the proportionality constant is a pure number of order 1 multiplied by d 2 . 
Therefore 

'C(X) ~ (33) 

We recall that the static potential is V — —h Q0 /2. The correction to the Schwarzschild metric 
is A e ff(x)h^ u (x), which computed at X gives a correction to the potential of the form 

V\X> ~ - GM ^'f (34) 
where A e ff < and "~" means that an adimensional constant of order 1 is omitted. 



4 The factor |A e// |<i 2 

Eq. (|34j) gives the first order correction to the 00 component of the Schwarzschild metric, due 
to a local cosmological term felt by the test particles. This correction is proportional to the 
product \A e ff\d 2 , where d is the size of the test particles. With the notation of Section 2, and 
supposing for the moment that the length scale £ which fixes the value of A e ff corresponds 
simply to d, the product \A e ff\d 2 can also be written as {d/p) 2 and is thus of order {l/d) a (eq. 

©)• 

The ratio (l/d) is very small, but since the exponent a is small, too, {l/d) a turns out 
to be a "reasonable" number also at a macroscopic scale; for instance, if d = 1 cm, then 
(l/d) a ~ (;lq-33^o.o3 ^ io _1 . Therefore the correction ([34]) is comparable in magnitude to the 



Newtonian component —GM/r (we set |X| = r in the following). In fact, the two contributions 
to /ioo are indistinguishable, since they have the same r-dependence, or in other words the 
cosmological correction just seems to lead to a renormalization of the factor GM. 



4.1 £ depends on the test particles' velocity 

The simple identification £ = d is not entirely correct, however. In the following, we show that 
£ has a weak dependence on the velocity of the test particle, and thus the r-dependence of the 
cosmological correction is not exactly 1/r. We first recall that the value of A e ff is the result of 
a spatial average. Now, suppose to follow the motion of the test body for a short time interval 
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At; the scale £ corresponds to the cubic root of the volume swept by the body in this interval, 
namely 

£~ (d 3 + d 2 vAt) 1/3 . (35) 

We see that £ depends also on the product vAt. For very large bodies (for instance, planets) 
one can consistently let At approach zero, so that the d 3 term in the bracket dominates, and 
£ = d with good approximation. For small test bodies, however, the term d 2 vAt will dominate, 
unless At is so small to correspond to an impracticable observation (for instance, if d = 0.1 m 
and v = 10 4 m/s, At should be less than 10~ 5 s). In other words, even though we are looking 
for the contribution of the local cosmological term to the potential V(r) at a given point, for 
small bodies A e ff depends on the volume swept in the motion during a short time interval and 
thus also on the velocity v. 



4.2 Dependence of V A on v(r) 

Let us go back to the expression (§) for A e //, namely |A e //| ~ (in this computation we 
can neglect a as compared to 2, thus 7 ~ 2). We set £ ~ (d 3 + d 2 v At) 1 ^ 3 like in eq. ([35]) and 
focus on the motion of a small test body, starting from a certain distance r . Disregarding 
the term d 3 in comparison to d 2 vAt we find that during this motion the adimensional ratio 
l^e//(^)|/|A e //(r )| is determined by the ratio of the velocities at the distances r and r : 



|Ae//(r)| _ (v{r) V 2/3 



(36) 



|A e //(r )| \v(r )J 

Note that At has dropped from this expression, as one should expect. Inserting it into eq. ([34] 
we obtain: 



V A (r) 



GM 



v(r) 



-2/3 



.™^H^J ' (37) 

According to the remarks at the beginning of this section, the constant factor d 2 \A e ff(r )\ 
is of order ~ 1. The r-dependence of the factor (v(r)/v (ro)) -2 ^ 3 causes a deviation from the 
pure 1/r behavior. This dependence can be obtained taking into account the total energy 
conservation. In the present approximation it suffices to consider the usual non relativistic 
conservation law: 

mMG 1 2/ . mMG 1 2/ . 

1 — mv (r) = 1 — mv (ro) (38) 

r 2 r 2 

or, denoting v 2 (r ) by Vq 

v 2 (r) 2MG f r Q -r \ . 2g Q r Q 

— — = 1 + — —[ =1 + — — r -r 39 

Vq Vq \ rr J Vq r 
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where go is the gravitational acceleration at the distance r . The constant go/v^ depends on the 
specific situation. For instance, the field of the sun at distance r = 1 AU is g ~ 1CT 3 m/s 2 , 
and taking as velocity of the test object vq ~ 10 4 m/s one obtains go/v" 2 ~ 10~ n m~~ l . 

For small displacements from ro, the general r-dependence in eq. (|39| ) can be expanded in 
powers of (r — r ) as follows (6 = go/v 2 ) 



, . , - . , -Kr-r ) + -6 26-- (r-r r + ... (40) 

After inserting this expression into eq. ([37])) one finds: (i) a potential 1/r, which corresponds as 
noticed above to an invisible "renormalization" of the product MG; (ii) an irrelevant constant 
potential; (iii) a potential linear in r, namely 

Vh n {r) ~ GM [rf 2 |A e// (r )|] b Ub - -) r. (41) 

Carrying on the expansion (|40|) to higher orders, one gets potentials growing like r 2 , r 3 ... This 
holds, however, only for r very close to r . In general, the expression (|39|) should be used. 



5 Discussion 

In the previous sections we investigated some consequences of the insertion into Einstein field 
equations of a local, scale-dependent cosmological term. By this we mean an effective cosmo- 
logical constant A e ff (essentially due to the quantum fluctuations of the gravitational field) 
whose strength depends on the size of the involved spatial region. 

For instance, if we consider the motion of a test particle in a given background field, the 
effective cosmological term felt by the particle depends on the volume swept by the particle in 
a short time interval during its motion. We computed the contribution of this term in the case 
of a Schwarzschild background, and found that the potential V = —h 00 /2 acting on a small test 
particle, usually of the Newtonian form V = —GM/r, is corrected by a small term depending 
on r as t> _2//3 (r)/r, being v(r) the particle velocity. 

The notion of a local cosmological term which "escorts" a particle in its motion is consistent 
with the definition of A e ff in lattice quantum gravity as the average curvature measured in a 
region of given size. However, the extrapolation of this notion from lattice theory and its 
implementation in the field equations could be properly justified only within a yet-to-come full 
quantum field theory of gravity. The situation reminds certain issues from the pre-QFT years 
early in this century, when quantum concepts were inserted by hand into classical field theory. 
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Our present approach allows, in some sense, to take into account the effect of the quantum 
gravitational fluctuations on a classical curved background. 

The conclusions of our computation, namely the existence of a tiny non-Newtonian compo- 
nent of the field depending on the size of the test particles, can be related to recent experimental 
results which called into play such a small "unmodeled" gravitational acceleration to explain 
the fine-detail motion of some space probes (while the planetary motion would not be affected). 

However our phenomenological approach does not allow to fix the numerical proportionality 
constants appearing in eq.s (|3"TD, (f4T|). Furthermore, the power expansion in eq. ( ffip is not suit- 
able for direct comparison with the experimental data. These originate from a complex fitting 
procedure, applied to the acceleration measured over long paths and based on conventional 
subtraction of the known parts of the field (with GM obviously constant). If the results of Ref. 
[0] will be confirmed, it may be worth to consider the inclusion of the potential (|37| ) in these 
fits. 

This work was partially supported by the A.S.P. - Associazione per lo sviluppo scientifico 
e tecnologico del Piemonte, Torino, Italy. 
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